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STUDYING ISM MAGNETIC FIELDS AND TURBULENT REGIMES FROM 

POLARIMETRIC MAPS 

D. Falceta-Gongalves, 1,2 A. Lazarian, 2 and G. Kowal 3 
RESUMEN 

Mapas polarimetricos han sido usados para caracterizar el campo magnetico en nubes moleculares. Sin embargo, 
cs diffcil determinar las propiedades tridimcnsionales de mapas de proyeccion de dichas regiones. Por esta razon, 
simulaciones numericas pueden ser usadas como pruebas de las mediciones polarimetricas, y eventualmente 
revelar mas acerca de la relation entre la turbulencia y las lmeas de campo magnetico. En este trabajo 
realizamos una serie de simulaciones MHD de nubes moleculares turbulentas y creamos mapas sinteticos de 
emision polarizada de polvo, en las cuales variamos la direction del observador. Determinamos la correlation 
de la intcnsidad de la emision y el grado de polarization para los modelos. Pudimos rcproducir el decaimiento 
del grado de polarization en zonas mas densas sin suposicion alguna de las propiedades dc la componente dc 
polvo. La anticorrclacion viene la simple cancelation de los vectores de polarization a lo largo de la visual. El 
efecto es amplificado dentro de zonas de alta densidad ya que la configuration del campo magnetico se torna 
mas compleja. Estudiamos la funcion de distribution dc probabilidad, el espectro de potencias, y funciones de 
estructura de los angulos de polarization. Estos analisis estadisticos revelan gandres diferencias dependiendo 
del regimen turbulento (es decir sub 6 super-sonico, sub 6 super- Alfvenico). Entonces, estos metodos pueden ser 
utilizados en observaciones polarimetricas para caracterizar la dinamica de las nubes moleculares. Presentamos 
tambien un metodo Chandrassekhar-Fermi modificado para obtener la magnitud del campo magnetico local. 
La formulation sugerida no muestra limitaciones respecto al regimen turbulento u orientation. 

ABSTRACT 

Polarimetric maps have been used for the characterization of the magnetic field in molecular clouds. However, it 
is difficult to determine the 3-dimensional properties of these regions from the projected maps. For that reason, 
numerical simulations can be used as benchmarks for polarimetric measurements, and eventually reveal more 
about the interplay of turbulence and the magnetic field lines. In this work we make a number of MHD numerical 
simulations of turbulent molecular clouds and created their synthetic dust emission polarization maps, varying 
the direction of the observer. We determined the correlation of intensity emitted and polarization degree for 
the simulated models. We were able to reproduce the decay of the polarization degree at denser regions without 
any assumption regarding the properties of the dusty component. The anti-correlation arises from the simple 
cancellation of the polarization vectors along the line of sight. This effect is amplified within denser regions as 
the magnetic field configuration becomes more complex. We studied the probability distribution function, the 
power spectrum, and the structure function of the polarization angles. This statistical analysis revealed strong 
differences depending on the turbulent regime (i.e. sub/supersonic and sub/super- Alfvenic). Therefore, these 
methods can be used on polarimetric observations to characterize the dynamics of molecular clouds. We also 
presented a modified Chandrashekhar-Fermi method to obtain the intensity of the local magnetic field. The 
proposed formulation showed no limitations regarding orientation or turbulent regime. 
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by large scale mag netic fields ( Schleuning 19981 : 



ICrutcher et al. 19991 ). However, it is still not com 



pletely clear what is the role of the magnetic field in 
the dynamics of the ISM and what is its effect on the 
star formation process. Also, the ratio of the mag- 
netic and turbulent energ y in these environments is 
a subject of controversy dPadoan fc Norlund 2002t 
Girart. Rao fc Marronei 2006 ). Polarimetric maps 



have been extensively used for the determination 
of the magnetic fields in several astrophysical envi- 
ronments. This technique represents the best tool 
to characterize the vector components of the mag- 
netic field parallel to plane of the sky. For galax- 
ies and the intracluster medium this technique can 
be used for polarizatio n of synchrotron emission 
(|Hildebrand et al.ll2000l ). 



For a given polarization map of an observed re- 
gion, the mean polarization angle indicates the ori- 
entation of the large scale magnetic field. On the 
other hand the polarization dispersion gives clues on 
the value of the turbulent energy. This, as a conse- 
quence, can be used to determine the magnetic field 
component in the plane of sky. Chandrasekhar & 
Fermi (1953) introduced a method (CF method here- 
after) for estimating the ISM magnetic fields based 
on the dispersions of the polarization angle and gas 
velocity. Simply, it is assumed that the magnetic 
field perturbations are Alfvenic and that the rms ve- 
locity is isotropic. 

A promising approach to test this method is to 
create two-dimensional (plane of sky) synthetic maps 
from numerically simulated cubes. Ostriker, Stone 
& Gammie (2001) performed 3D-MHD simulations, 
with 256 3 resolution, in order to obtain polarization 
maps and study the validity of the CF method on 
the estimation of the magnetic field component along 
the plane of sky. They showed that the CF method 
gives reasonable results for highly magnetized me- 
dia, in which the dispersion of the polarization angle 
is < 25°. However, they did not present any other 
statistical analysis or predictions that could be use- 
ful to determine the ISM magnetic field from obser- 
vations. Polarization maps from numerical simula- 
tions can also be used in the study of the correlation 
between the polarization degree and the total emis- 
sion intensity (or dust column density). Observa- 
tionally, the polarization degree in dense molecular 
clouds decreases with the total intensity as P cx I~ a , 



with a = 0.5 — 1.2 (jGoncalves. Galli fc Walmslev 
20051 ). Padoan et al. (2001) studied the role of 
turbulent cells in the P vs I relation using super- 
sonic and super- Alfvenic self-gravitating MHD sim- 
ulations. They found a decrease of the polariza- 
tion degree with total dust emission within grav- 
itational cores, in agreement with observations, if 
grains are assumed to be unaligned for Ay > 3. 
When the alignment was assumed to be independent 
on Ay, the anti-correlation was not observed. Re- 
cently, Pelkoncn, Juvela & Padoan (2007) extended 
this work and refined the calculation of polarization 
degree introducing the radiative transfer properly. In 
that work, the decrease in the alignment efficiency 
arises without any ad hoc assumption. The align- 
ment efficiency decreases as the radiative torques be- 
come less important in the denser regions. However, 
it is still not clear the role of the magnetic field topol- 
ogy and the presence of multiple cores intercepted by 
the line of sight on the decrease of polarization de- 
gree. 

In this work we attempt to extend the previously 
cited studies improving and applying the CF method 
to different situations. For that, we used both sub 
and super-Alfvenic models, to study the role of the 
magnetic field topology in the observed polarization 
maps. 

2. NUMERICAL SIMULATIONS 

The simulations were performed solving the set 
of ideal MHD equations, in conservative form, as fol- 
lows: 
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dpv 
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(1) 



= pf, (2) 



— - V x (v x B) = 0, 



(3) 



with V ■ B = 0, where p, v and p are the plasma 
density, velocity and pressure, respectively, B is the 
magnetic field and f represents the external accelera- 
tion source, responsible for the turbulence injection. 
For molecular clouds, we may assume that the ratio 
of dynamical to radiative timescales is very large. 
Under this assumption, the set of equations is closed 
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with an isothermal equation of state p = c 2 p, where 
c s is the speed of sound. The equations are solved 
using a second-order-accurate and non-oscillatory 
scheme, with periodic boundaries, as described in 
Kowal, Lazarian & Beresniak (2007). 

Initially, we set the intensity of the x-directed 
magnetic field B ex t and the gas thermal pressure 
p. This allows us to obtain sub-Alfvcnic or supcr- 
Alfvenic, and subsonic or supersonic models. 

The turbulent energy is injected using a random 
solenoidal function for f in Fourier space. This, in 
order to minimize the influence of the forcing in the 
formation of density structures. We inject energy at 
scales k oc L/l < 4, where L is the box size and / is 
the eddy size of the injection scale. The rms velocity 
SV is kept close to unity, therefore v and the Alfven 
speed va = B/^/AjTp will be measured in terms of 
the rms SV. Also, the time t is measured in terms 
of the dynamical timescale of the largest turbulent 
eddy (~ L/SV). 

We performed four computationally extensive 3D 
MHD simulations, using high resolution (512 3 ), for 
different initial conditions, as shown in Table 1. We 
simulated the clouds up to t m ax ~ 5, i.e. 5 times 
longer than the dynamical timescale, to ensure a full 
development of the turbulent cascade. We obtained 
one subsonic and three supersonic models. One of 
the supersonic models is also super- Alfvenic. Each 
data cube contains information about parameter- 
ized density, velocity and magnetic field. As noted 
from Eqs. (1) and (2), the simulations are non self- 
gravitating and, for this reason, the results are scale- 
independent. 

Regarding the gas distribution in each model we 
found an increasing contrast as we go to higher sonic 
Mach number, independently on the Alfvenic Mach 
number. Subsonic turbulence show a gaussian distri- 
bution of densities, while the increased number and 
strength of shocks in supersonic cases create smaller 
and denser structures. In these cases, the density 
contrast may be increased by a factor of 100 - 10000 
compared to the subsonic case. The magnetic field 
topology, on the other hand, depends on the Alfvenic 
Mach number. Sub- Alfvenic models show a strong 
uniformity of the field lines, while the super- Alfvenic 
case shows a very complex structure. Both effects, 
the density contrast and the magnetic field topology, 
may play a role on the polarimetric maps, as shown 



TABLE 1 

DESCRIPTION OF THE SIMULATIONS 



Model 


P 


Bext 


Description 


1 


1.00 


1.00 


subsonic & sub-Alfvcnic 


2 


0.10 


1.00 


supersonic & sub- Alfvenic 


3 


0.01 


1.00 


supersonic & sub- Alfvenic 


4 


0.01 


0.10 


supersonic & super- Alfvenic 



further in the paper. 

3. RESULTS 

3.1. Polarization Maps 

To create the polarization maps we assumed the 
dust polarization to be completely efficient and that 
the radiation is originated exclusively by thermal 
emission from perfectly aligned grains. Under these 
assumptions, the local angle of alignment (tp) is de- 
termined by the local magnetic field projected into 
the plane of sky, and the linear polarization Stokes 
parameters Q and U are given by: 



q = p cos 2tp sin 2 i, 
u = p sin 2ip sin 2 i, 



(4) 



where p is the local density and i is the inclina- 
tion of the local magnetic field and the line of sight. 
We then obtain the integrated Q and U, as well 
as the column density, along the LOS. Notice that, 
for the given equations the total intensity (Stokes 
/) is assumed to be simply proportional to the col- 
umn density. The polarization degree is calculated 
from P = \/Q 2 + U 2 /I and the polarization angle 
4> = a tssx{U/Q). 

In Figs. 1 and 2 (see Falceta-Gongalves, Lazarian 
& Kowal [2008], FLK08 hereafter), we show the ob- 
tained maps of column density and the polarization 
vectors for Models 3 and 4, respectively, for a line 
of sight perpendicular to the original magnetic field 
orientation. In Fig. 1, the polarization vectors are 
mostly uniform, and parallel to the external mag- 
netic field. Here, fluctuations on the polarization 
angle are seen within the condensations, where the 
kinetic pressure overcomes the magnetic pressure. In 
Fig. 2, we show the column density and polarization 
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Fig. 1. Polarization of emission and column density maps for Model 3 (Ms ~ 7.0 and Ma ~ 0.7) with H ex t perpendicular 
to the line of sight. The complete map (512x512 pixels) (Upper-left) and the zoomed regions (100x100 pixels). The 
sensitivity in simulated observations is assumed to be 0.3 of the maximum emission. Here, regions where the signal is 
less than 0.3 do not show polarization vectors, and P m ax = 97%. 
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Fig. 2. Polarization of emission and column density maps for Model 4 (Ms ~ 7.0 and Ma ~ 7.0) with H ex t perpendicular 
to the line of sight. The complete map (512x512 pixels) (Upper-left) and the zoomed regions (100x100 pixels). The 
sensitivity in simulated observations is assumed to be 0.3 of the maximum emission. Here, regions where the signal is 
less than 0.3 do not show polarization vectors, and P m ax = 97%. 
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maps of the super- Alfvenic case (Model 4). Here, 
the kinetic energy is larger than the magnetic pres- 
sure almost everywhere. As a consequence the gas 
easily tangles the magnetic field lines. The angu- 
lar dispersion is larger and the polarization degree 
is smaller when compared to the sub-Alfvcnic case. 
For the super-Alfvenic case, the orientation of the 
magnetic field regarding the LOS is irrelevant to the 
polarization maps. The strongly magnetized cases 
(Models 1, 2 and 3) present very similar polariza- 
tion maps. Also, we must note that for LOS parallel 
to the original magnetic field, the polarization maps 
seem similar to the super-Alfvenic case. 

Regarding the polarization degree, obser- 
vations have shown that a decorrelation be- 
tween the polarization degree and the col- 
umn densities is detected for several objects 
dMatthews fc Wilson! l2002t lLai. Girart fc Crutcherl 
120031 ; IWolf. Launhardt fc Hennindl2003h . Typical 
the polarization degre e follows the relation P oc J~ 7 , 
wher e 7 ~ 0.5 — 1.2 I Goncalves. Galli fc Walmslevl 
2005J) and I is the total intensity. Usually, it is 



assumed that the grain alignment efficiency may 
be the general cause of this effect. However, as 
noticeable from Figs. 1 and 2, the polarization 
degree is smaller within high column density regions 
for all models. Note that we assumed perfect 
grain alignment, independent on the local density. 
Therefore, the decorrelation in our simulations must 
have a different cause. Another possibility could be 
the effect of the averaging along the LOS. Denser 
regions present a less uniform magnetic field. The 
annihilation of polarization vectors perpendicular 
to each other could cause the decrease in the 
polarization degree. 

In Fig. 3 we show the correlation between the 
polarization degree and the column density for the 
different models, assuming that the magnetic field is 
in the plane of the sky. For all models, the polariza- 
tion of high column densities tend to decrease to the 
minimum value (~ 20%P max ). This minimum polar- 
ization degree should be zero for homogeneous den- 
sity and random magnetic field. In inhomogencous 
media it depends on the number of dense structures 
intercepted by the line of sight. The major contri- 
bution for the polarized emission comes from dense 
clumps, which only are few along the LOS. This poor 
statistics results in a non-zero polarization degree. 



o. 
o. 




Fig. 3. Correlation between averaged polarization degree 
and the column density for the different models (FLK08). 
P m ax is 100%, 98%, 97% and 85% for Model 1, 2, 3 and 
4, respectively (FLK08). 

We also noticed that in subsonic turbulence the po- 
larization degree is large even for the higher column 
densities. It occurs because in the subsonic models 
the contrast in density is small and the simulated 
domain is more homogeneous. Also, the number of 
dense clumps, which are able to tangle the field lines, 
is reduced in the sub-sonic case. On the contrary, 
for the super-Alfvenic case we obtain the minimum 
polarization degree everywhere, i.e. purely random 
magnetic field components. For the decreasing part 
of the plots, we obtained a best fit with a correla- 
tion exponent 7 = 0.5. If grain alignment was prop- 
erly implemented in the calculations, the value of 7 
should increase. 

3.2. Statistics of Polarization Angles 

The histograms of polarization angles are shown 
in Fig. 4. In the upper panel we show the histograms 
for the sub- Alfvenic (Models 1, 2 and 3) and the 
super-Alfvenic (Model 4) cases, with the mean mag- 
netic field lines perpendicular to the LOS. The polar- 
ization angles present very similar distributions and 
almost equal dispersion for the sub-Alfvenic cases. 
This happens mainly because they do not depend on 
the density structures, but on the magnetic topol- 
ogy. Strongly magnetized turbulence creates more 
filamentary and smoother density structures (i.e. low 
density contrast) if compared to weakly magnetized 
models and, most importantly the magnetic field 
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Fig. 4. Histograms of polarization angle of the different 
models with 9 = (up), and for Model 3 and different 
magnetic field orientations in respect to the line of sight 
(angles 6) (bottom) (FLK08). 



lines are not highly perturbed. For Model 4, the dis- 
tribution is practically homogeneous, which means 
that the polarization is randomly oriented in the 
plane of sky. It occurs because the turbulent /kinetic 
pressure is dominant and the gas is able to easily 
distort the magnetic field lines. 

In the bottom panel of Fig. 4 we show the polar- 
ization angle histograms obtained for Model 3 but 
for different orientations of the magnetic field. The 
dispersion of the polarization angle is very similar for 
inclination angles 8 < 60°, and increases for larger 
inclinations. It may be understood if it is noted that 
the projected magnetic field B s k y = B ext cos 9 is of 
order of the random component SB. It shows that 
the dominant parameter that differ the distributions 
of <fi is the uniform magnetic field projected in the 
plane of sky, and not the intensity of the global mag- 
netic field. 

Another statistical analysis, which helps us to 



understand the topology of the polarization angles, 
is based on structure functions. The second order 
structure function (SF) of the polarization angles is 
defined as the average of the squared difference be- 
tween the polarization angle measured at 2 points 
separated by a distance I: 



(5) 



SF(Z) = (|0(r + l)-4>(r) 



The structure functions calculated for the differ- 
ent models are shown in Fig. 5. In the upper panel we 
present the SFs obtained for Model 3 with different 
values of 0. As expected, all curves present a pos- 
itive slope showing the increase in the difference of 
polarization angle for distant points. However, the 
small scales part of the SF presents a plateau ex- 
tending up to I r~j 4 — 5 pix. This range corresponds 
to the dissipation region and may also be related to 
the smallest turbulent cells. Surprisingly, the log-log 
slopes arc equal independently on the inclination of 
the LOS (for 1 between 3 and 20 grid points). The 
same behavior was obtained for all models. In the 
bottom panel we show the SFs calculated for the 
different models with 8 = 0. It is noticeable the in- 
crease in the SF for higher Mach numbers. However, 
the slopes are notably different. The maximum slope 
is a ~ 1.1, 0.8, 0.5 and 0.3 for Models 1, 2, 3 and 
4, respectively. Observationally, the molecular clou d 
M17 shows a ~ 0.5 up to I = 3pc (|Dotson|[l996h . 
which would be in agreement with a cloud excited 
by supersonic and sub (or critically)- Alfvenic turbu- 
lence. 

3.3. The Chandrasekhar- Fermi Method 

Chandrasekhar & Fermi(1953) proposed a 
method to estimate the ISM magnetic fields based 
on the dispersion of polarization angles and the rms 
velocity. Basically, it is assumed that the magnetic 
field perturbations are Alfvenic, i.e. Sv cx SB^fp. If 
the rms is isotropic, we may relate both as: 



£U L os 



(6) 



If one wants to expand the applicability of the 
CF method for cases where the random component 
of the magnetic field is comparable to the uniform 
component, or for larger inclination angles, it is nec- 
essary to take into account two corrections. Firstly, 
we must introduce the total magnetic field projected 
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Fig. 5. Structure functions of polarization angle for 
Model 3 with different magnetic field orientations regard- 
ing the line of sight (angles 8) (up), and for the different 
models with magnetic field perpendicular to the line of 
sight (6 = 0) (bottom) (FLK08). 



in the plane of sky £? sky ~ + SB, where B^ 
represents the mean field component projected on 
the plane of sky. We assume here, for the sake of 
simplicity, that SB is isotropic. Wc assume that 
the SB/B is a global relation and, in this case, we 
may firstly obtain the dispersion of <j) and then cal- 
culate its tangent. Substituting S<fi in Eq. (6) by 
tan((5</>) ~ SB/Bgky, we obtain the modified CF 
equation: 

BX+SB^V^P^y (7) 

Since the CF-method strongly depends on the 
resolution (bacause of the dispersion of angles and 
velocities), we applied Eq. (7) to our simulated 
clouds, taking into account the effects of finite res- 
olution. We calculated the average of the density 
weighted rms velocity along the LOS (SV\ os ) and the 
dispersion of the polarization angle (S<p) within re- 
gions of R x R pixels. To simulate a real cloud we 
substitute the parametric values of the Model 3 for 
?i H = 10 3 cm.- 3 and T = fOK, and B ~ 50/iG. The 
results showed a correlation between the obtained 
magnetic field from Eq. (7) and the resolution, 
which may be described by: 

B cf = B° CF (l + ^ ) , (8) 

where R represents the observational resolution (to- 
tal number of pixels), C and B^ F are constants ob- 
tained from the best fit. Bq F represents the value of 
Bcf for infinite resolution observations, i.e. the best 
magnetic field estimation from the CF method. 

The fit parameters, as well as the expected val- 
ues of the magnetic field from the simulations for all 
models, are shown in Table 2. Here, the magnetic 
fields are given in units of the mean field B cxt ■ Since 
the simulations are scale independent, one could 
choose values of B cxt to represent a real cloud, in 
accordance with the parameters of Table f . As an 
example, assuming a cloud with njj = 1 3 cm -3 and 
T = fOK, and = O.Of (Model 3), we get B ~ 50/xG. 
Choosing differently the density, temperature or the 
model given by the simulations, i.e. (3, we obtain a 
different mean magnetic field. The obtained param- 
eter C is very similar for the different inclinations, 
but are different depending on the model, mainly be- 
cause it is related to the scale on which the dispersion 
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of the polarization angle changes. Since C seems to 
depend on the model and not on the inclination it 
could also be used by observers to infer the physical 
properties of clouds from polarization maps. 

As a practical use, observers could obtain polari- 
metric maps of a given region of the sky for different 
observational resolutions (e.g. changing the resolu- 
tion via spatial averaging). Using the CF technique 
for each resolution and, then apply Eq. (8) to de- 
termine the asymptotic value of the magnetic field 
projected into the plane of sky Bq F . 

4. CONCLUSIONS 

In this work we presented turbulent 3-D high res- 
olution MHD numerical simulations in order to study 
the polarized emission of dust particles in molecular 
clouds. We obtained synthetic dust emission polar- 
ization maps calculating the Stokes parameters Q, U 
and I assuming a perfect grain alignment and that 
the dust optical properties are the same at all cells. 
Under these conditions, we were able to study the 
polarization angle distributions and the polarization 
degree for the different models and for different in- 
clinations of the magnetic field regarding the LOS. 

As main results, we obtained an anti-correlation 
between the polarization degree and the column den- 
sity, which is in agreement with observations, with 
exponent 7 ~ —0.5. This value is related to random 
cancellation of polarization vectors integrated along 
the LOS, while larger indices require extra physics, 
e.g. dependency of the dust alignment efficiency with 
the local density. 

We showed that the overall properties of the po- 
larization maps are related to the Alfvenic Mach 
number and not to the magnetic to gas pressure 
ratio. Also, we studied the PDF's and structure 
functions of the polarization angles, which showed 
a degeneracy of the results with the Alfvenic Mach 
number and the angle between the magnetic field 
and the LOS. Zceman measurements of dense clouds 
may be useful to help remove this degeneracy as it 
could provide the magnetic field component along 
the LOS and, therefore, the inclination. 

Finally, we presented a generalization of the CF 
method, which was showed to be useful for: (i) the 
determination of the total magnetic field projected 
in the plane of sky, and (ii) the separation of the two 
components -B s ky and SB; 
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FALCETA-GONgALVES, LAZARIAN, & KOWAL 



TABLE 2 
CF METHOD ESTIMATES 



Model 


en 


C 




DCXt / JD 

^sky/^cxt 


Btot/B ext 


3 





20 ±5 


1.24 ±0.09 


1.00 


1.25 


3 


30 


24 ±5 


0.98 ±0.08 


0.87 


1.11 


3 


45 


25 ±5 


0.78 ±0.07 


0.71 


0.96 


3 


GO 


33 ±5 


0.48 ±0.05 


0.50 


0.75 


3 


90 


31 ±5 


0.26 ±0.03 


0.00 


0.24 



1 7±5 0.97 ±0.08 1.00 1.11 

2 10 ±5 1.07 ±0.07 1.00 1.16 
4 34 ±5 1.18 ±0.07 1.00 1.41 



